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We consider single-inclusive forward pion production in high-energy proton-proton collisions at 
RHIC energies. A good baseline description of the transverse momentum distributions at high 
rapidity is obtained within Mueller's dipole formalism with an anomalous dimension incorporating 
an "extended geometric scaling" window between the saturation and DGLAP regimes. We then 
formulate pion production for transversely polarized projectiles within the same approach. We 
assume that an azimuthal, spin-dependent asymmetry arises from the so-called Sivers effect and 
■ investigate the single transverse-spin asymmetry An at 200 GeV and 500 GeV center-of-mass energy. 

A simple parameterization of the Sivers functions from the literature compares reasonably well with 
the high-energy STAR data if the overall normalization is scaled up by at least a factor of two. The 
STAR data might therefore indicate that the Sivers effect is significantly stronger than thought so 
far. We also analyze higher- twist contributions to An and find that they largely cancel. 
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Large single transverse-spin asymmetries (SSAs) in p^ p — > it X have been observed in fixed target experiments 
over the last 15 years. A non-vanishing SSA implies that the azimuthal distribution of pions in the final state 
depends on the direction of the transverse spin of the polarized projectile proton. These observations pose a 
challenge for theorists who want to explain such asymmetries in terms of quark and gluon degrees of freedom. 
It is clear that the asymmetries result from spin-orbit couplings and chiral dynamics in low-energy QCD, but to 
describe such couplings and chiral symmetry breaking within quantum field theory is far from straightforward. 
Attempts to describe the observable 



^ da(p1 p — > 7T X) — da(p^ p — > it X) 
dcr(pT p — » 7T X) + dcr(p-L p — > 7r X) 

perturbatively at leading twist and within the framework of collinear factorization where partonic transverse 
momenta in hadrons are assumed to play a subdominant role (because the associated scale is of order Aqcd), 
result in very small asymmetries 0. Several proposals that go beyond leading twist or collinear factorization 
have been put forward 0. Q|. 0. Q and can in principle account for the large asymmetries. They all involve new 
polarization-dependent parton distribution or fragmentation functions, which however cannot be calculated within 
perturbative QCD (pQCD) because of the soft dynamics. So far, these functions could also not be calculated 
from Euclidean QCD on a lattice. Hence, reliable predictions can only be made once the unknown functions have 
been extracted from experiments. As a further complication the evolution properties of the functions, needed to 
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describe the energy dependence of the asymmetry observables, are also not known yet, because of the much more 
complicated operator structure as compared to the ordinary parton densities. Despite these limitations, a fairly 
successful phenomenology has emerged. Given the large magnitude of the effects it is certainly worth exploring 
such SSAs further, both experimentally and theoretically. 

Due to the rather low energies and transverse momenta probed by fixed target experiments, a factorized pQCD 
description of the process may be somewhat questionable. Therefore, the high-energy collider experiments with 
polarized protons at BNL's Relativistic Heavy-Ion Collider (RHIC) were eagerly awaited. Recently, the STAR 
collaboration at RHIC has found that the large single spin asymmetries persist at collider energies of y/s — 200 GeV 
and transverse momenta up to pr ~ 2.5 GeV/c 7]. This has been confirmed by the BRAHMS experiment at 
RHIC ; 8j. Ref. Q suggested that the asymmetries at low and high center of mass energies are actually two 
different phenomena, which was based on the observation that only at RHIC energies the unpolarized inclusive 
pion production cross section is described satisfactorily by a leading- twist, next-to-leading order (NLO) pQCD 
calculation. The fact that this does not work for fixed-target energies casts some doubt on the validity of the 
factorization assumption for those cases. 

In this paper we shall specifically address the high-energy collider data on SSAs. As a first step, we attempt to 
obtain a good baseline for inclusive pion production in unpolarized pp collisions. SSAs are usually measured at 
forward rapidities y (in the hemisphere of the polarized projectile) where they are found to be largest. The large 
rapidity coverage of the forward tt° detectors at STAR enable measurements of pion production up to y ~ 4 t 10j . 
At large y, however, the kinematics is very asymmetric: the hadrons in the final state emerge from collisions of 
projectile partons with large light-cone momentum fraction x\ ~ 0.1 - 1 with target partons carrying a very small 
momentum fraction x% -C 1. When a s log l/x2 is of order 1, small- a; effects in the target become important and 
modify its gluon distribution function as compared to that obtained in the DGLAP approximation (specifically, 
its anomalous dimension). The purpose of this paper is to show how such small- x effects can be accounted for 
in the description of SSAs at forward rapidities and high energies. Their relevance for deep inelastic scattering 
(DIS) at HERA and forwar d pion p roduction in pA collisions at RHIC has already been discussed in (among 
others) refs. GHEEl and [Il|lEllllI3> respectively. 

We will assume that the asymmetry arises predominantly from the so-called Sivers effect Q. Within the 
framework of spin and transverse momentum dependent distribution and fragmentation functions, it has been 
shown 0| that other possible mechanisms, such as the Collins effect cannot describe the asymmetries without 
other additional contributions, most importantly the Sivers effect. We therefore simply restrict to the Sivers effect, 
but inclusion of the other types of spin effects will proceed along similar lines. 

The first study of small- a; effects, more specifically of saturation effects within the so-called "Color Glass Con- 
densate" (CGC) formalism, on a sing le spin asymmetry concerned production of polarized A baryons in collisions 
of unpolarized protons on nuclei |19| ■ No effects from small- a; quantum evolution, such as geometric scaling vio- 
lations (see below), were taken into account though. Also, we shall show here that higher- twist saturation effects 
play no significant role for pp collisions at RHIC energy, as may be expected from estimates of the saturation 
scale Q s . 

This paper is organized as follows. In the next section we describe a model for the Sivers effect using a well-known 
parameterization taken from the literature. In sec. Ill- A we formulate single-inclusive forward hadron production 
in unpolarized pp collisions within the CGC formalism and show that good agreement with the experimental data 
emerges. Next, by applying the formalism to p^p collisions, in sec. III-B we compare our results for the SSA to 
recent STAR data at large rapidity (~ 3.8). Sec. IV is devoted to a discussion of the importance of higher-twist 
effects and finally a summary is given in sec. V. 

II. A MODEL FOR THE SIVERS EFFECT 

The Sivers effect is a correlation between the direction of the transverse spin of the proton and the transverse 
momenta of its unpolarized partons. The Sivers effect in the process p^ p — » ir X has first been analyzed by 
Anselmino et al. |5j. The Sivers functions for valence u and d quarks were extracted from a fit to fixed-target 
data under the assumption that the asymmetry arises solely from this effect. That fit relied on a calculation of 
both the polarized and unpolarized cross sections within the LO DGLAP pQCD approach and partial inclusion 
of partonic transverse momentum dependence. In the present paper we will analyze whether the Sivers function 
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from jf| is consistent with the high-energy STAR data, when incorporated in the CGC formalism. To be precise, 
we shall use the Sivers functions presented in refs. [ill EI- Although updated fits [23 incorporate more detailed 
transverse momentum dependence, the simpler parameterizations of [2(J,|21| will suffice for our present purposes. 
The number densities of q = u, u, d, d, ■ ■ ■ partons in a polarized proton are defined as 

fq/pi (X, h) = f q /p(x, kt) + ~A N f q/p1 (x, k t ) S * ' ( 2 ) 

1 \St\\P\\k t \ 

where St and P are the transverse polarization and three-momentum vectors of the proton, respectively, and k t 
is the quark's intrinsic transverse momentum. The function f q / p {x,kt) is the density of quarks with transverse 
momentum k t in an unpolarized proton, from which the familiar quark distribution f q / p (x) arises by integration 
over kt, see below. A N f q / p i (x,k t ) denotes the Sivers function (see ref. for a comparison to another frequently 
used notation f^j, 0). 




FIG. 1: Illustration of the p^ p — » ir° X process. The kinematics is simplified by selecting the pion production plane (x — z) 
perpendicular to the polarization St = (0, S y , 0) of the proton. Parton transverse momenta are taken to be (anti-) parallel 
to the x-direction. We define left (right)-hand production of pions in the direction of St X P > (< 0) to match the 
experimental situation. 



Fig. ^ illustrates the p tt° X collision process, simplifying the kinematics by taking the tt° production 
plane (x — z) perpendicular to the polarization of the proton, St = (0, S y , 0), where the Sivers effect is maximal. 
The collision axis is along the z-direction, i.e. P — (0,0, P z ). Following the standard experimental definition, we 
define the left and right directions of ir° production in the forward (positive hadron rapidity y^) region from the 
perspective of the polarized proton as indicated in the figure. This definition implies that if A N f q / p r {x, k t ) > 0, 

then the quark in the polarized proton has a preference for moving to the left, i.e. with kt pointing in the positive 
x-direction. 

The difference between upward and downward spin polarizations of the proton, which enters in the numerator 
of An, therefore picks up a leading contribution from the Sivers functions. Flipping the spin of the proton is 
equivalent to changing the sign of k t : 

AN fq/ P r (x, k t ) = f q/p i (x, k t ) - f q/p i (x, k t ) = / g / p T (x, k t ) - f q/p r (x, -kt) 

= A f q/p ,{x,k t ) , (3) 

I <s* I \P\ \h\ 

where A N f q / p i (x, kt) is an odd function of k t . The Sivers effects cancel in the sum over polarization states, which 
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enters in the denominator of A jy and yields (twice) the leading-twist unpolarized function, 

f q / p r{x,kt) + f q /pi(x,k t ) = 2f q/p (x,k t ). (4) 

Note that / d 2 k t f q / p (x, k t ) = f q / p (x). 

The inclusion of intrinsic fct-dependence in the parton distributions is not straightforward from a theoretical 
point of view. First of all, factorization theorems involving transverse momentum dependence are generally harder 
to derive. Considerable progress has been made for specific processes ;25J which are, however, less complicated 
than the one considered here. We shall view the factorized formulas employed below as reasonable phenomeno- 
logical extensions of collinear factorization at the semi-hard observed transverse momenta we are considering. 
Furthermore, we refer to ref. for a derivation of collinear factorization from a one-loop analysis of unpolarized 
parton scattering off a CGC target. 

Another important theoretical issue regarding Sivers functions is the calculable process dependence |2fj. Since 
the fitted functions from refs. apply to the same process as considered here (and to the same partonic 

subprocess in the sense that qg scattering is dominant in the considered kinematic region), we will not worry 
about this issue here. 

As already mentioned above, the Sivers functions appear in the difference between upward and downward spin 
polarizations of the proton 

da(p ] p~> hX) - dcr(p l p^ hX) cx J d 2 k t [/ g/pT (x, k t ) - f q/pl (x, k t )\ eg) da qp ^ q ' x (k t ) ® D h/q , (z h ) 

= J d 2 k t A N f q/p ,(x,kt)®da qp ^ qX (k t )®D h/q ,(z h ) , (5) 

where da qp ^ q x (kt) denotes the cross section of a parton scattering off an unpolarized proton. The fragmentation 
functions Dh/ q i(zh) describe the fragmentation of a final-state parton to a hadron with momentum fraction 
Zh- Note that the Sivers functions correspond to distributions of unpolarized partons and so the elementary 
hard-scattering cross sections are the usual expressions for unpolarized partons, in both the numerator and the 
denominator of An . 

Because the convolution integrand has to be an even function of kt, while the Sivers function is odd, one finds 
J d 2 k t A N f q/p r (x, k t ) <g> da qp ^ q ' x (k t ) 
= ( I + f ) d 2 k t A N f q/p r(x,kt)®da q P- q ' x (k t ) 

\Jxk t >0 Jx-k t <Cl/ 

= f d 2 ktA N f q/p ,(x,kt)®[do qp ^ qX (kt)-do qp ^ qX (-k t )] 1 (6) 

Jh.p. 

where "h.p." indicates that the integration extends over only the positive half plane with x ■ k t > 0. 

Following joi 1 2 ("| we assume that the Sivers function A N f q / p i (x, k t ) is sharply peaked about an average trans- 
verse momentum = \J (k 2 ) such that 

A N f q/p , (x, k t ) ~ A N f q/p i (x) S^(k t - k%) . (7) 

The unit vector k t is taken to lie in the positive half of the scattering plane, x ■ k t > 0. Then 

J d 2 k t A N f q/p ,(x,k t )®da q ^ qX (k t )®D h/ql (z h ) 

~ A N f q/p , (x) ® [da q P^ q ' x (kfkt) - da q P^ q ' x (-k° t k t )} ® D h/q ,(z h ). (8) 

Before the STAR data was published, the only available quantitative knowledge of the Sivers functions 
A N f q /pi (x) was due to fits to SSA data from fixed-target pp experiments and to low energy semi-inclusive 
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DIS data, cf. ref. |27| and references therein. Although we will briefly comment on some more recent fits below, 
for this paper we will mainly restrict ourselves to the simpler parameterizations that resulted from fixed-target pp 
data. We only allow for a different overall normalization of the Sivers functions. For valence quarks, the functions 
were parameterized as follows poll^ 1 : 

A N f q/pl (x) = K Af MM Nq x ; (i _ x) h . (9) 

The x-dependent average transverse momentum of valence quarks in the polarized proton is taken as 

k°(x) = 0.47 x - 68 (1 - x) 0A8 M p , (10) 

for all flavors. M p = 1 GeV denotes the proton mass. The maximal intrinsic transverse momentum of rj 210 MeV 
is on the order of the inverse proton radius and therefore reasonable. At large x, one does not expect much higher 
intrinsic transverse momenta. The remaining parameters are given by 

N u<d /A = (3.68, -1.24), a u , d = (1.34,0.76), b u>d = (3.58, 4.14) . (11) 

It was pointed out in ref. \2\\ that the original values of N q from pol ] were misquoted and have to be multiplied 
by a factor of 4. In addition, in © we allow for further adjustment of the magnitude of the Sivers functions by 
a factor f if required by the new high-energy data from STAR. 



III. FORWARD PARTICLE PRODUCTION AT HIGH ENERGIES 



A. Unpolarized pp collisions 

The forward pion spectra in pp collisions obtained by the STAR Collaboration at RHIC (y/s = 200 GeV) [2jJ 
are generally known to agree reasonably well with leading-twist NLO pQCD calculations in the DGLAP ap- 
proximation |9j, |3Jj . Although this approach is recognized as providing a much better description of semi- hard 
particle production at RHIC than at lower energy [[j , not all commonly employed sets of parton distribution and 
fragmentation functions (FFs) perform equally well. To be specific, the NLO pQCD result with CTEQ6M PDFs 
and KKP FFs is consistent with the data for rapidity (r]) ~ 3.8, but compares less well at (77) = 4.0, except in 
the high-p t region. The Kretzer FFs lead to better agreement with lower-p t data, presumably due to its smaller 
gluon fragmentation to pions as compared to the KKP set, which may not be realistic though and one needs more 
precise extractions of the pion FFs at moderately high-p t (pt = 1 ~ 3 GeV/c) before definite conclusions can be 
drawn. The fact that some standard PDFs and FFs fail to give a good description at forward rapidities may find 
a natural explanation in the observation that these data are sensitive to small- a; effects in the target. 

Hence, we would like to analyze these high-energy, forward processes at RHIC within another theoretical 
approach, namely the "Color Glass Condensate" (CGC) formalism. Contrary to the DGLAP approach, it resums 
small- a; effects in the gluon distribution function. The Sivers functions, on the other hand, are unaffected as the 
SSA arises from large- x valence quarks on the projectile side. 

The CGC approach is most commonly applied to the study of gluon saturation effects at sufficiently small x 
and/or Q 2 |l5l l31|. At a given x, the transverse momentum below which the unintegrated gluon distribution 
saturates due to non-linear terms in the QCD evolution equations is called the saturation momentum Q R (x). It 
grows approximately as a power of energy NLO BFKL evolution gives Q s (x) ~ l/x x with A ~ 0.3 34]. 

HERA phenomenology indicates that Q s ~ 1 GeV at x ~ 10 -4 

The deep saturation (non-linear) regime is therefore of minor relevance to semi-hard {pt ^ 1 GeV) particle 
production in pp collisions at RHIC energy. This will be verified in a more quantitative way in section IIVI 



1 We note that these expressions are understood to apply at a fixed (average) scale. The scale dependence is not known yet and will 
not be considered here. 



6 



However, the CGC formalism is also applicable to the so-called "extended" geometric scaling (EGS) regime |32T . l35| 
which emerges above the saturation line up to transverse momenta of about Q gs = Q 2 /A 3> Q s , where A is a 
nonperturbative soft scale expected to be of order Aqcd- In deep inelastic scattering (DIS) geometric scaling 
implies that at low x (<, 0.01) the proton structure function depends only on the scaling variable r = Q 2 /Q 2 s (x) 
rather than on Q 2 and x separately. 

Fits to DIS data from HERA flsTl found that approximate geometric scaling holds not only at small Q 2 , 
where this is expected due to saturation effects, but over a much broader window Q 2 < Q 2 < Q 2 S , thereby 
providing evidence for the existence of such a new regime. It arises if the solution of the LO-BFKL evolution 
equation is expanded to second order around the saturation saddle point The first term gives geometric 

scaling while the second "diffusion" term contributes to scaling violations 2 . Due to this diffusion term, the 
anomalous dimension 7 governing the small- a; evolution of the gluon distribution with rapidity y = log(l/x) is 
shifted from 7 S ~ 0.63 (the BFKL saddle point in the vicinity of the saturation line) by A7 cx \og(l /r t Q s )/y, 
where r t ~ 1/Q denotes the dipole size. As Q approaches ~ Q gs from below, violations of geometric scaling grow 
to order unity and the anomalous dimension reaches 7dglap ~ 1 (more precisely, 1 — 7dglap = 0(a s )). Ref. 0] 
finds that the average shift required for a good fit to the HERA data is not small, (A7) ~ 0.2. Thus, the behavior 
of the anomalous dimension which determines the (geometric) scaling violations plays a key role for high-energy 
ep scattering. 
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FIG. 2: "Phase diagram" of the target proton at RHIC energy (*J~s = 200 GeV). yh and pt are the rapidity and the transverse 
momentum of the produced hadron, respectively. The dashed line depicts the boundary Q 3 (yh) of the saturation region; 
the dotted line that of extended geometric scaling, Q gs {yh)- For this plot we fixed the projectile parton's momentum 
fraction to the typical value x\ — 0.1. The solid line corresponds to the boundary of phase space, xf = 1. 



The same approach has also been applied to high-p t ( >, 1 GeV) hadron production in d Au collisions at RHIC 
over a wide interval of rapidity yh flT) ] . Indeed, the data was shown to agree rather well with a parameterization 
of the anomalous dimension which increases logarithmically with p t from 7 = "f s to its asymptotic value 7 ~ 1, 
while decreasing with y as A7 ~ l/y at very large rapidity. The d Au data from RHIC is therefore consistent with 



2 To find Q gs ~Q 2 /A one should in fact estimate Q gs from the transition point between the LLA and DLA saddle points, respectively, 
rather than from the diffusion term in the expansion of LO-BFKL about the saturation saddle point. For details see, for example, 
section 2.4.3 in ref. and references therein. 
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the presence of an extended geometric scaling regime above the saturation region which widens towards forward 
rapidities (i.e. smaller x 2 )- The EGS window is characterized by the need to resum logarithms of \jx% rather than 
of Q 2 . However, the twist expansion remains valid, all- twist resummation being necessary only in the saturation 
regime Q 2 < Q 2 . 

We illustrate the various kinematic regimes emerging in pp collisions at RHIC in Fig. The boundaries of 
these regions were determined as follows. The saturation regime of the target proton extends to 



(x 2 ) = lGeV 



(s) 



A/2 



(12) 



where, as already referred to above, xo = 3.0 x 10~ 4 and A = 0.3 llj. X2 is related to the rapidity yh of the 
produced hadron by kinematics: x\jx 2 = exp(2j//,) (see, for example, appendix B in For this figure, we 

fixed X\ to a value typical for a valence quark, x\ = 10 _1 . The EGS regime corresponds to the window between 
Q s and Q gs , where the anomalous dimension 7 (given below in eq. l(THj| ) grows from "f s to ~ 1. In fig. [3 its upper 
limit is taken to be Q gs — Q 2 /A with A = 0.5 GeV, for illustration. Above Q gs , 7 ~ 7dglap and so the DGLAP 
approximation applies; note that this of course includes the DLL regime of large yh and large pt- The boundary 
of phase space, finally, is reached when the Feynman-a; of the produced hadron xf = (pt/y/s) exp(y? l ) = 1. 

This figure confirms that the all-twists saturation regime plays no prominent role at RHIC energy and pt ^ 1 
- 1.5 GeV/c, regardless of rapidity. However, for yh^3 the EGS window opens wide and eventually covers 
essentially the entire semi-hard regime of particle production as the DGLAP region is cut off by energy-momentum 
conservation constraints. This implies that at such rapidities, 7 does not reach 1 at any accessible transverse 
momentum. It may therefore be plausible that small- a; effects are relevant for semi-hard particle production in 
the forward region of pp collisions at RHIC. 

Forward particle production leads to a strong asymmetry in the momentum fractions x\ ~ 1 and 12 < 1 which 
are probed in the projectile and target, respectively. The process can then be described using an asymmetric 
DGLAP (g) CGC factorization scheme developed in ref. [l6[ . 



dN h 



axi — 



fq/p(xi,Pt)N F (~Pt' X ^) D h/q 



dyhd 2 pt (2"7r) 2 J Xf x f 

+ fg/p(xi,Pt)N A ( — p t ,x 2 ) D h/g (—,Pt 



(13) 



The indices q in f q / p and D h / q are implicitly summed over all participating quark flavors. In what follows, the 
mass of the produced hadrons is neglected because -C pt, which allows us to equate to a good approximation 
the pseudorapidity r\ with the rapidity yh- 

The approach l|13|l includes DGLAP evolution of the projectile parton distribution functions f q / p (xi ) p 2 ) and of 
the fragmentation functions D h / q (xF /xi,p 2 ) with the factorization scale Q 2 = p\ . It resums small- a; logarithms 
in the target as well as DGLAP logarithms of Q 2 on the projectile side and has been derived to leading order in 
a s |lfij |. Exact rather than small- a; approximated splitting functions are to be employed, such that the projectile 
momentum is conserved, which is important at large x\. On the other hand, the target proton is described by the 
dipole forward scattering amplitudes Np,Ai where Np corresponds to a projectile quark impinging on the target's 
small- a; gluons, while Na applies to a projectile gluon. These arise |3jj as averages of two-point functions of 
Wilson lines in the corresponding representation (running up and down the light-cone, separated by a transverse 
distance r y ) o ver the gluon field of the hadron [13. In the high-energy limit one therefore recovers Mueller's dipole 
approach 38]. For a recent comparison of this approach with HERA diffractive electroproduction data we refer 
to ref. 

The dipole profile employed here is similar to that from the fit presented in ref. jljj to HERA DIS data, except 
for the presence of the anomalous dimension 7 < 1 [lj, [Tj, [l]| [l(| [l7| : 



N A (r t ,x 2 ) = 1 - exp 



^2)) 



(14) 
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rt denotes the transverse size of the dipole. Np is obtained here by the replacement Q 2 S — > (Cp /Ca) Q 2 = 
(4/9) Q 2 , corresponding to a weaker coupling of a quark relative to a gluon projectile to the target field. Due to 
this color factor, our effective saturation momentum entering Np is somewhat smaller than the original Golcc- 
Biernat-Wiistoff (GB-W) fit [HI but agrees with previous comparisons to RHIC data [bH llii Il7| and with the 
analysis of HERA DIS data presented in ref. 0] . 

The Ansatz 1|14J) exhibits saturation at r t > 1/Q S and also reproduces the Na ~ f\ DGLAP limit as rt — ► 
(assuming that also 7 — > 1 in that limit) 3 . It also displays exact geometric scaling when 7 is constant since then 
N(rt, x) depends only on the scaling variable p = r t Q s - Violations of geometric scaling arise due to the dependence 
of 7 on rt, increasing from 7 S at r t ~ 1/Q S to 7dglap ~ 1 at very short distances. Here, the anomalous dimension 
7 (ft j £2) of the gluon distribution is parameterized as 



i(n,x 2 ) = 7s + (1 - is) 



\og(l/r*Ql(x 2 )) 



At/2 + + log(l/r t 2 Q2( X2 )) 



(15) 



with 7 S ~ 0.627; 1/2 = log \ jx% is (minus) the rapidity of the target partons. This function agrees with the 
solution of LO-BFKL evolution with saturation boundary conditions in the saddle point approximation [13( , and 
with the general theoretical limits mentioned above. However, the parameterization (|15(l also includes sublcading 
corrections ~ d^Jy^ which govern geometric scaling violations at subasymptotic rapidities. With d ~ 1.2 the 
anomalous dimension l|15|) was shown to provide a good description of dAu data from RHIC over a wide range 
of rapidity We note that 1)15(1 applies only in the EGS regime, while 7 = j s stays fixed within the saturation 
region (which is, however, not important here, see section ITvj) . 
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FIG. 3: Transverse momentum distributions of forward inclusive ty° from unpolarized pp collisions at y/s = 200 GeV. The 
STAR data [29| were taken at average pseudorapidities of (77) = 3.8 and 4.0, respectively. Theoretical results correspond to 
two different sets of parton distribution and fragmentation functions, i.e. CTEQ5-LO+KKP-LO (I) and MRS-G+BKK-LO 
(II), respectively, all evaluated at Q 2 — p\. They give x 2 P er data point of 0.9 (0.3) at y = 3.8 and 0.7 (1.0) at y = 4.0 for 
set I (II). 



3 Note that in order to reproduce the ~ l/<7t DGLAP power-law tail in momentum space one needs to introduce an additional 
logarithm of 1/Vj in the exponent of 1141 . For 7 = 1, NA(lt) then drops approximately exponentially in the high-density regime 
about qt ~ Qs, t urni ng into ~ LO DGLAP behavior as qt — * 00 lilt |4(1 . A similar modification of the GB-W model was 

considered in ref. l4ll . 
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To determine the transverse momentum distribution of hadrons via eq. (|13|l requires taking the Fourier trans- 
forms of the dipole profiles l|14|l to momentum space. To simplify this task, we shall replace 7(7*4, X2) by 7(1/54, X2), 
where qt — (x\/x F )pt is the Fourier conjugate of r t . This provides a good approximation at large y^ (small X2) 
since there 7 increases rather slowly with q t . 

In fig. we present our numerical results for the ^-distributions of forward inclusive pions from pp collisions 
at = 200 GeV. They correspond to rapidities yh — 3.8 and 4.0, respectively. We employ two different sets 
of parton distribution (PDF) and fragmentation (FF) functions: CTEQ5-LO PDFs H2 and KKP-LO FFs 0, 
which is nowadays a common combination. The other set is MRS-G PDFs 0] and BKK-LO FFs as used 
in ref. [2(il | to extract the parameterization of the Sivers functions from experimental data obtained at lower 
energies. We observe rather good agreement of both sets in the interval 1.2<,p t <, 2.4 GeV/c with recent STAR 
data from RHIC [29| . For the overall normalization, we need to multiply our LO results by a factor K ~ 1.4 for 
the CTEQ5-LO+KKP-LO set I, while K ~ 1.0 for the MRS-G+BKK-LO set II. The X 2 per data point is 0.9 
(0.3) at y = 3.8 and 0.7 (1.0) at y — 4.0 for set I (II), which is rather good. Since both sets work equally well 
within error bars, we will henceforth employ set II which does not require a K-factor. 

The results from fig. [3] correspond to a scale Q 2 = p\ . The scale dependence in eg. ()13f) arises entirely from the 
LO PDFs and FFs and amounts to a weakly p t -dependent normalization factor [0|. If the scale dependence of 
the LO Sivers functions turns out to be similar to that of the unpolarized PDFs, such an overall JC-factor would 
cancel in the ratio of the polarized to unpolarized processes. The evolution properties of the Sivers functions are, 
however, largely unknown at present and as a consequence we have to postpone a more detailed study of such 
issues to the future. 



B. p^p collisions and single transverse-spin asymmetry 

The single-inclusive pion distribution of eq. I|13|) is written in a collinearly factorized convolution form, da(pp — > 
hX) <~ fqi p <g> N F <S> D h / q . We shall formulate the single transverse-spin asymmetry An in an analogous way but 
partly include transverse momentum dependence in order to incorporate the Sivers effect. For the numerator of 
eq. this results in 

do{p"p -> hX) - da(p l p -> hX) cx J d 2 k t {f q / p i (xi, k t ) - f q / p i{x 1 , k t )] ® N F (x 1 ,q t ,k t ) eg) D h / q (x F /xi) 

~ A N f q/p i( Xl ) ® [N F (x u q t - k°) - N F (x u q t + k°)}® D h/q (x F / Xl ) ,(16) 

where we used eq. © to arrive at the last line. Also, we simplified the kinematics (as was done in refs. j2£j and 
therefore has to be considered as part of the fit of the Sivers function) by picking up only contributions where k t 
is either parallel or anti-parallel to q t . 

On the other hand, the denominator of is given by 



da(p ] p —> hX) + da{p l p — > hX) cx J d k t [f q/p r (xi, k t ) + f q / p i (xi, k t )} ® N F (x\, q t , k t ) eg) D h/q {x F /xi) 

~ 2f g/p (x 1 )®N F (x 1 ,q t )®D h/q (x F /x 1 ), (17) 

where eq. (@J and the fact that f q / p {x,kt) is an even function of k t have been used. Thus, when displaying 
explicitly all arguments of the various functions the asymmetry becomes 



A I \ val-g 



E A^/pt^) g) [N F (^p t - kly 2 ) - N F (f^ + kly 2 )] ®D h/q 



2 E f q / P (xi,P 2 t) ® N F (jfcpt, y 2 ) ® D h/q (|f + f g / P (xi,P 2 ) ® N A (Urft.lte) ® D h/g (|f ' 

(18) 

Only valence-quark contributions are accounted for in the numerator, according to the parameterization of the 
Sivers functions from ref. pfij . On the other hand, the unpolarized cross section in the denominator includes 
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the contributions of all active quark flavors and of gluons. At large rapidities, however, it is also dominated by 
valence quarks, ft is evident from ljl8|l that a good description of the elementary quark q t distribution and of its 
derivative is required to reliably extract the Sivers function from An- 
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FIG. 4: Single transverse-spin asymmetry An in the rapidity interval yh = 3.8 - 4.2 for yfs = 200, 500 GeV, with -Ka/ = 2 
(two times larger Sivers functions than ref. |H). Th e STAR data were taken at an average pseudorapidity of (rj) = 3.8 
and y/s = 200 GeV. 



In Fig. 01 we compare our numerical results for the asymmetry An to STAR data Q taken at y/s = 200 GeV. 
We also show predictions for larger rapidities, y% — 4.0 and 4.2, and for higher energy, s/s = 500 GeV. The 
theoretical curves were obtained from eq. I|18|l with the MRS-G [44| PDFs for the unpolarized projectile proton 
and the BKK-LO FFs, the dipole l|14|) . the anomalous dimension (|15|l . and the saturation momentum (I12|l . 
For the Sivers functions, we adopted the parameterization I|9I10I11|I with -K"a/ = 2, corresponding to a two times 
larger Sivers effect as compared to the original fit from ref. 20] . 

One observes that our curve at yh — 3.8 goes up more slowly than the data and peaks at large pt around 
Pt ~ 3.5 GeV/c. At higher rapidity, however, the curves increase more steeply and the peak shifts to smaller 
p t . At fixed rapidity, the asymmetry decreases rather rapidly with energy. This is mainly a kinematic effect 
in that the typical momentum fraction x\ of projectile partons decreases, and so does the Sivers function 
Given the theoretical uncertainty in the Sivers functions and the sizable errors of the data, the agreement can be 
considered reasonable. Even with the present large error bars, however, the STAR data appear to imply larger 
Sivers functions than previously thought |2fJ (as reflected by K^f — 2), at least if An indeed arises entirely from 
the Sivers effect. If more precise data becomes available in the future, our formalism (extended to allow for more 
general k t dependence of the Sivers functions if needed) could be used for a complete refit of the Sivers functions. 

A few additional comments regarding the factor K&f ^ 2 found here are in order. The parameterization (|9I11|) 
of the Sivers function resulted from a fit to low-energy fixed-target data, for which an accurate description of the 
pion transverse momentum distribution (and of its slope) is hard to obtain. It is therefore not entirely surprising 
that the new high-energy data, while confirming qualitative features of the parameterizations, may still require 
some quantitative adjustments. Of course, such adjustments would alter predictions for asymmetries of other 
processes, like for instance the Drell-Yan process [22|. When high energy data with smaller error bars becomes 
available we intend to perform a more detailed investigation of the x- and fcf-dependence of the Sivers functions, 
rather than of their overall normalization only. 

We should also point out that a good description of the STAR data can also be achieved by increasing the 
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intrinsic transverse momentum kt (x) of valence quarks substantially, as recent parameterizations show |46j . This 
is because the Sivers distribution is proportional to the intrinsic transverse momentum, and another factor 
of k t arises in eq. ((161) from the difference of the contributions of quarks with k t parallel or anti-parallel to the 
momentum transfer; see for example eq. (|23() below. However, the process considered here is dominated by large 
x on the projectile side (x\ ~ 0.1 - 1), where intrinsic transverse momenta should probably not exceed the inverse 
proton radius by much. This is consistent also with our finding (and that from NLO DGLAP approaches) that 
the forward inclusive pion distribution from unpolarized collisions can be described rather well without assuming 
large intrinsic transverse momenta of projectile valence quarks. 



IV. HIGHER-TWIST EFFECTS ON SSA 



The unitarized dipole profiles Np,A entering 1)16117(1 resum all higher-twist effects which arise at low transverse 
momentum. It is interesting to disentangle the leading-twist contribution to the observables considered here. For 
Pt > Qsi one expects that higher-twist contributions are suppressed. 

We estimate the leading contribution by introducing 



iV>(n,j/2;c) 



1 — exp 



(ArtQi(y 2 )/9) 



■y(r t ,y 2 ) 



(19) 



For c = 1 we recover the full Np(rt, y 2 ) from (f 1 4|) . In the limit c — > 0, however, the Fourier transform of 1(19(1 
isolates the leading contribution in the expansion with respect to Q s /qt- For q t ^ one obtains 



N F (qt ^0,y 2 ;c) 



2- 



d 2 r t e^N F (r t ,y 2 ;c) 

dr t rtMnqt) exp [-£(4r t 2 Q2 (2/2 )/ 9 )7(n,y 2 ) 



(20) 



We know analytic forms of the Fourier transform 1(20(1 only for 7=1/2 and 1 [l^. In the limit c — > one finds 



Nl- 1/2 {q u y 2 -c) 



Np (q t ,y 2 ;c) = 



c^q vr ( Q 



(Q s c/6) 2 + q 2 



Z ' 2 3Q2 V qt 



9?r 

3!? exp 



9 f qt 



4c \Q 



(21) 
(22) 



Hence, for 7 = 1/2, representing the LO-BFKL anomalous dimension without saturation boundary conditions, 
the Fourier transformed dipole profile exhibits a leading-twist tail ~ (Qs/qt) 3 (we recall that ~ (Qs/qt) 4 emerges 
in the DGLAP regime). On the other hand, the GB-W model 11] corresponding to 7 = 1 shows no power-law 
tail at large q t . It arises entirely from all- twist resummation and has no expansion in powers of Q s /qt- 

For the anomalous dimension from 1(151) the Fourier transform and its c — ► limit cannot be obtained analytically. 
We therefore show numerical results in Fig. OD Here, we set yh — 3.8, pt — 2 GeV/c, which represent typical 
values, and take x\ = xf such that Pt = qt- 

One observes that the profile for 7 = 1/2 is almost constant in the interval < c < 1 and approaches a finite 
value as c — > 0, which agrees with the analytic result ((21(1 . This confirms that for these kinematic conditions the 
dipole with anomalous dimension from LO-BFKL without saturation boundary condition is dominated entirely 
by the leading-twist contribution. Next, we consider the q t -dependent anomalous dimension ((15(1 which generates 
violations of geometric scaling and interpolates to the DGLAP regime at high q t . For this case, the Fourier 
transformed dipole profile decreases slowly as c — ► from its value at c = 1. This decrease corresponds to higher- 
twist contributions of about 8%. Thus, higher- twist contributions to single-inclusive forward hadron production 
in pp collisions are not very large once an anomalous dimension with the described features is taken into account. 
Quantitative comparisons to the anomalous dimension from leading-twist NLO DGLAP approaches for the present 
kinematics would be interesting. 
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FIG. 6: Same as Fig.^but now for the difference N F (q t - k°) - N F (q t + k°). 



Next, wc address the same question for the difference appearing in the numerator of the asymmetry An which 
essentially corresponds to the derivative of Np(q t ). For this quantity, the leading contribution is suppressed by 
an additional power of k® /qt as can be easily realized by considering the leading contribution for 7 = 1/2 as an 
example: 



3 1.0 

6^A + . 



(23) 
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Due to this additional factor, the magnitude of Np(q t — h®) — Np(q t + fc°) as shown in fig- El is smaller than that 
of Np(q t ) from fig. El Here, too, the c — > limit isolates the leading power of Q s /qt from this quantity. The curve 
for 7 = 1/2 is again flat over the entire interval, implying leading-twist dominance. On the other hand, for 7 = 1 
(GB-W) the leading-twist contribution vanishes completely. For the (^-dependent anomalous dimension l|15fl we 
observe higher- twist contributions of about 11%. 
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FIG. 7: The influence of the anomalous dimension on the asymmetry Am at yh = 3.8 and 1 < p t < 3 GeV/c. 

Finally, Fig. [7| depicts the asymmetry An over the range 1 < pt < 3 GeV/c at = 3.8 and y/s — 200 GeV 
for 7 = 1/2, 1, and 7 = 7dhj from eq. (|15(1 . For the latter case, we also show A/v as obtained from the leading- 
twist contribution to Np alone. Since the full and leading-twist curves are almost identical, we conclude that 
higher- twist contributions in the numerator and denominator cancel to a large extent. On the other hand, it 
turns out that the growth of the anomalous dimension with transverse momentum (geometric scaling violations) 
does give rise to a slightly steeper A^ipt) as compared to a constant 7 = 1/2. The GB-W dipole, which does 
not permit a twist expansion, gives a very small asymmetry. The reason for this behavior is the unrealistically 
large contribution from gluons to the unpolarized cross section in the denominator of An at high rapidity: in 
momentum space, NA(qt) ~ exp[— q 2 /Q 2 ] falls off less rapidly with q t than Np(q t ) ~ exp[— (Ca/Cf) Qt/Q'sl- 

V. SUMMARY 

In this paper, we have analyzed recent high-energy data from STAR on inclusive pion production at forward 
rapidities and on single transverse spin asymmetries (SSA) in p^p collisions. Particle production at large rapidities 
involves very small momentum fractions in the target, hence we propose and use a factorization approach which 
accounts for small- a; effects on the gluon distribution of the target. Within this "Color Glass Condensate" 
approach, the target is described by unitarized two-point functions of light-like Wilson lines (dipole forward 
scattering amplitudes) rather than the usual leading- twist, leading-log (DGLAP) gluon distribution function. 

We have shown explicitly that the deeply non-linear ("saturation") regime is of little relevance for particle 
production in pp collisions at RHIC energy. The dipole scattering amplitude is therefore well approximated by its 
leading-twist limit N(r t ,x) ~ {r t Q s {x)) 2l ^ rt ' x \ Nevertheless, small-x evolution not only leads to saturation for 
transverse momenta q t < Q s (x) but also to a so-called "extended geometric scaling" window above Q Sl which has 
been confirmed by fits to HERA DIS data. Within the EGS window extending from Q s {x) to Q gs (x) ~ Q 2 s (x)/ \ 3> 
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Q s (x), the anomalous dimension 7 of the gluon distribution function increases smoothly from its value at the 
BFKL saddle-point with saturation boundary conditions, 7 S ~ 0.63, to its value in the dilute DGLAP regime, 
7dglap ~ 1- Once this feature is incorporated, a very good baseline for forward pion production in unpolarized 
pp collisions is obtained (x 2 per data point as 1). A comparison of our 7 to that from NLO DGLAP approaches 
in the kinematic regime relevant for STAR would be interesting, given that the predicted pion ^-distributions 
are similar. 

Building on this baseline, we then proceed to analyze single transverse-spin asymmetries in the forward region. 
We restrict ourselves to SSAs originating from the Sivers effect, which is nowadays believed to give the dominant 
contribution within the framework of transverse spin and transverse momentum dependent parton distribution 
functions. The approximately 10% contributions from higher twists to the polarized and unpolarized cross sections, 
respectively, largely cancel in their ratio An, which is therefore very well approximated by the leading- twist 
contribution. Furthermore, we find that a parameterization of the Sivers functions obtained previously from a fit 
to fixed-target pp data j2lj within the LO DGLAP approach provides a reasonable description of the high-energy 
STAR data, if their overall normalization is scaled up by at least a factor of two. This might indicate that the 
Sivers effect is significantly stronger than thought so far and that small- a; effects do not cancel in the ratio, by 
affecting the slope of the cross section. More quantitative statements about the Sivers functions must await data 
with smaller errors and over some range of energies. 
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